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Abstract. A system of two species of fermions of different mass confined in a one¬ 
dimensional harmonic trap is studied with an exact diagonalisation approach. It is 
shown that a mass difference between fermionic species induces a separation in the 
density of the lighter flavour independently of the number of particles. The mechanism 
behind the emergent separation is explained phenomenologically and confirmed by 
direct studies of the ground state of the system. Finally, it is shown that the separation 
driven by a mass difference, in contrast to the separation induced by a difference of 
populations, is robust to the interactions with thermal environment. 


1. Introduction 

Development of cooling and trapping techniques of a few quantum particles (bosons as 
well as fermions) have opened completely a new area of modern quantum engineering 
[1-5]. In todays experiments, it is possible to rigorously control not only the strength of 
mutual interactions between atoms but also the precise number of quantum particles of 
different flavours conhned in a trap of almost arbitrary shape [5-7]. These possibilities 
are strongly related to the fundamental theoretical questions about the properties of a 
mesoscopic number of particles which up to now had not been answered in a satisfactory 
way. One of these questions is about the transition from two-, three-body physics [8-10] 
that are quite well understood with standard approaches, to macroscopic systems, that 
are well described with the help of sophisticated methods of condensed matter physics. 
A hrst attempt to give an experimental answer to this question was given in [11] where 
the emergence of the Fermi sea with increasing number of fermions was studied in a 
few-body system. In the attractive regime the effective appearance of two-body pairing 
as a function of the system’s size was observed [12] and a theoretical description of this 
phenomenon was given [13-16]. Recently, many different results in the area of ultra-cold 
few-body systems were reported both in experimental [17-24] and theoretical [25-47] 
works. 
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In a recent paper [48] , a theoretical analysis of phase separation in a one-dimensional 
imbalanced mixture of two kinds of fermions was made. It was shown that strong 
interactions with a single impurity dipped in a mesoscopic Fermi sea lead to a separation 
of density in the bulk flavour. This separation is directly induced by the imbalance in the 
number of fermions of the two types. In particular it can not be observed in a balanced 
system. However, it is known that in balanced systems consisting of flavours of different 
atomic mass, separation is also possible. For example, in the recent paper [49], density 
separation driven by a mass difference was shown in a system conhned in a uniform trap. 
Since in experiments with ultra-cold fermionic mixtures it is possible to conhne atoms 
with different masses [6,7,21,23,32], therefore it is interesting to study the properties of 
such a system in detail. In this paper, we study the separation mechanism in a mixture 
of two-flavour fermionic atoms of different mass conhned in a one-dimensional harmonic 
trap. The whole analysis is performed with an exact diagonalisation of the many-body 
Hamiltonian. 

At this point it is worth noting that Hubbard-like models with imbalanced mass of 
particles were also intensively studied [50,51]. 


2. The Model 


In this paper we study a two-havour system of a few ultra-cold fermions conhned in a 
one-dimensional harmonic trap of frequency u and interacting via a short-range delta¬ 
like potential. The fermions of diherent havours are fundamentally distinguishable and 
in principle may have diherent masses ^ m^. If one expresses all quantities in the 
natural units of the harmonic oscillator related to one of the havours (in the following we 
choose a =4,) then the Hamiltonian of the system in the second-quantized form reads: 

«= 5 : 




dx T|(a:)T|(a:)T 4 .(a;)T^(a:) 


( 1 ) 


where To-(x) annihilates a fermion with spin cr at a point x and the dimensionless 
interaction constant g = {mi/fiidujY^'^giY) is related to the ehective one-dimensional 
interaction gi-Q between fermions of opposite spins. Experimentally this ehective one¬ 
dimensional model is obtained by applying very strong harmonic conhnement in two 
remaining perpendicular directions [17]. Then it can be assumed that the dynamics 
in these directions is frozen and particles occupy only the ground state of the trap. 
A direct relation between ehective interaction giD and the s-wave scattering length 
which can be tuned with Feshbach resonance was discussed with all the details in [28]. 
The single-particle Hamiltonians take into account the mass ratio a = m^/mi and 
have the following forms 

1 d^ a 9 


1 d^ 1 2 
“ 2d^ ^ 2^^ ’ 


It is worth noticing that in principle both havours of fermions can be conhned in traps 
with diherent frequencies. Here we assume that the experimental parameters are tuned 
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such that the frequencies are equal. At the level of the single-particle Hamiltonian the 
only difference between flavours is related to the mass difference. This assumption has 
some consequences on the properties of the system. For example, the single-particle 
excitation energy is flavour-independent and it is equal to huj. 

Since particles of different flavours are distinguishable, the held operators 
for opposite havours fulhl standard commutation relations T'f(x), 


T-t-(a;), = 0. However, the held operators for the same havour obey the 

canonical fermionic anti-commutation relations ~ ^') 

f = 0. Obviously, the Hamiltonian (1) commutes separately with the 


number of fermions of a given havour = fdx Therefore, further analysis 

will always be performed under the assumption that the number of fermions in a given 
havour is known and well dehned experimentally. This assumption is reasonable since 
in ongoing experiments with few fermion systems the number of particles is controlled 
to a very high accuracy [5,11]. 


3. Method 

The properties of the system described by the Hamiltonian (1) are studied with the 
exact diagonalisation approach what was previously adopted for equal mass fermionic 
mixtures [14,15,52], in both the repulsive and attractive case. First, we decompose 
the held operators To-(a;) onto the basis of the eigenfunctions of the appropriate single¬ 
particle Hamiltonian, i.e., il/o-(x) = J2i4>ii7{x)dia. Note, that due to the mass diherence, 
the basis functions of diherent havours have diherent spatial spread and they are related 
by the scaling relation 0q(x) = (l)ii{y/a x). In this language, the Hamiltonian (1) 
has the following form 

n = EE ^icF “t“ ^ ^ ^ijkl (3) 

(7 i ijkl 

where hjo- = The single-particle energy Ei = {i + 1/2) and the two-body 

interaction term Uijki can be calculated directly 

Uijki = g j dx(t)*^{x)(t)*^{x)ct)ki{x)(t)i^{x). (4) 

Since the wave functions (j)ia{x) form a complete set in the space of single-particle states, 
the Hamiltonian (3) is fully equivalent to the original Hamiltonian (1). In practice, we 
have to cut the decomposition of the held operators at some sufficiently large level 
index Amax- For appropriately chosen cut-oh, the results do not depend on its value. 
For example, for the parameters = 3, = 1, a = 1 and 5 ^ = 10 it is sufficient to 

take Amax = 18 single-particle eigenstates in the decomposition of the held operators. 
Then, the dimension of the many-body Hilbert space is dim{H) = 14688. 

After obtaining the matrix form of the Hamiltonian (1) in the Fock basis for a given 
number of fermions and interaction strength g, and after subsequent diagonalisation, we 
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Figure 1. (a) Spectrum of the Hamiltonian (1) for = 3 and = 1 as a function 
of the interaction g. Depending on the ratio of atomic masses a the spectra differ 
qualitatively for strong repulsion. For example, the degeneracy manifolds in the limit 
of g —> oo have different dimensions, (b) Spectrum of the Hamiltonian (1) for = 3 
and = 1 and strong interaction 5 = 10 as a function of the mass ratio a. It is seen 
that for increasing a, the energy of the two many-body eigenstates belonging to the 
ground-manifold rapidly split and the ground-manifold has only a two-fold degeneracy. 


find the M lowest many-body eigenstates |Gj) and corresponding eigenenergies £i of the 
Hamiltonian. 

4. Spectrum of the Hamiltonian 

The first significant consequence of the mass difference between flavours is visible in the 
spectrum of the Hamiltonian. In the standard situation, when a = 1, the Hamiltonian 
of the system is symmetric under exchange of the flavours. As a matter of fact, this 
symmetry, together with the Pauli exclusion principle, leads directly to the appearance 
of a degeneracy in the manifolds of many-body eigenstates of the Hamiltonian in the 
limit of infinite repulsion [52,53]. Each degenerate manifold is spanned by states with 
the same spatial distribution of particles but with different characters of the irreducible 
representation of the exchange-symmetry group mentioned above. For large but finite 
repulsions the spectrum of the Hamiltonian is not degenerate but it factorises into distant 
manifolds of quasi-degenerate eigenstates. The mechanism that lifts the degeneracy and 
its experimental consequences have been broadly discussed recently [39,52], 

It is worth noticing that the case of = 1 is solvable analytically for any 

value of the mass ratio a. The two-body problem is brought to a single-particle one by 
introducing centre of mass and relative distance coordinates. In these coordinates, 
the Hamiltonian of the system is decomposed into two single-particle problems for 
which analytical expressions for eigenenergies and eigenfunctions are known [54], Many 
properties of two particles conhned in harmonic trap were recently studied in fermionic 
as well as in bosonic systems [17,55-59]. 

Whenever a ^ 1, the exchange-symmetry is violated since the interaction part of 
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the many-body Hamiltonian, in contrast to the single-particle one, is not symmetric 
under exchange of particles. This directly leads to different spectral properties of the 
Hamiltonian in the limit of large repulsion. To illustrate this mechanism, in Fig. la we 
compare spectra of the Hamiltonian (1) with a = 1 and a = 40/6 (corresponding to the 
mass ratio in the mixture of lithium and potassium atoms [6,7,21,23,32]) calculated for 
a hxed number of fermions = 3 and = 1 as a function of interaction g. As it can 
be seen, for large interactions the degeneracy manifolds do split. Our numerical results 
suggest that this mechanism is present for different and N^. Unfortunately, due to 
the numerical complexity, our method enables us to study only those conhgurations for 
which < 10. However, for all of them the same mechanism is present if the 

mass ratio a is large enough. Therefore, we believe that the mechanism is very general 
and it is present for any number of particles. 

It is quite natural that for different pairs of the degeneracy may be lifted into 
sub-manifolds with different numbers of eigenstates. To give a better explanation of 
how this degeneracy is lifted we show in Fig. lb the spectrum of the Hamiltonian (1) 
for = 10 as a function of the mass ratio a. For a = 1 the ground-manifold is spanned 
by four many-body eigenstates. However, for a larger a the eigenenergy of the two 
states rapidly grows. Due to changes of the shapes of the corresponding single-particle 
wavefunctions the contact interaction between flavours becomes signihcant. It is 
seen that, the same mechanism of quasi-degeneracy splitting repeats also for higher 
manifolds. 

5. Properties of the ground-state 

The mass difference between fermionic species signihcantly changes not only the 
properties of the spectrum of the Hamiltonian but also some properties of many-body 
states. First, let us discuss changes of the many-body ground-state |Go). For a = 1 the 
properties of the ground-state were studied recently in [48]. It was shown that repulsive 
interactions between species lead to phase separation whenever the number of particles 
is imbalanced, i.e. ^ N^. In particular, the authors studied the interaction of a 
single impurity dipped in the Fermi sea of the remaining species > 1 and = 1). 
They found that for strong interaction, the bulk flavour is pushed out of the centre of 
the trap and two spin domains precursors appear. The remaining impurity is located in 
the middle with an almost unchanged density prohle. This observation can be visualized 
by plotting the single-particle densities dehned as follows 

p^{x) = (Go I I Go). (5) 

Note, that the density is normalized to N^j. In the upper panel of Fig. 2 we show 
density prohles for = 3 and = 1 obtained with the exact diagonalisation approach 
for different g. In this way we reconstruct the corresponding result of [48] by using an 
exact diagonalisation method which is accurate near g = ^ for weak interaction. 

The situation is qualitatively different when a difference of mass between flavours 
is introduced. As an illustration, we repeat all the calculations for a = 40/6 (lower 
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Figure 2. Single-particle densities (thin blue line, heavier flavour) and Pi{x) 

(thick red line, lighter flavour) calculated in the ground-state of the system for different 
interactions g and mass ratios a = 3 and = 1). For weak interactions, the 
single-particle densities reproduce the density profile for ideal gasses in both species. 
When interaction is strong enough, a phase separation takes place in the system. For 
equal masses, the separation occurs in the species with the larger number of atoms 
(as predicted in [48]). On the contrary, for a = 40/6 the lighter species undergoes 
separation. We have checked that this behaviour is independent on the number of 
fermions in individual flavours (see main text and Fig. 3). 


panel in Fig. 2). Obviously, in the limit of vanishing interaction the only difference 
comes from the fact that the heavier impurity is better localised in the middle of the 
trap. However, for strong enough interactions the separation is not present in the bulk 
subsystem, but a signihcant separation of the lighter impurity appears. Note that in this 
case, the separation of the lighter flavour means that the impurity is in the ” Schrodinger 
cat”-like state, i.e. the probability of hnding the particle is non-zero in the two, spatially 
separated regions of space. It was reported recently that this sort of split state can also 
be seen for impurities in a bosonic bath [60]. Also a quite similar effect of the phase 
separation driven by the mass difference was studied recently for spin-1/2 fermions in a 
uniform trap [49] . 

The mechanism of separation induced by the mass imbalance is fundamentally 
different form separation driven by interactions. For example, for equal masses a = 1, 
the separation is present in the ground-state only for different number of particles 

7 ^ N^) and it appears always in the flavour with larger N„. It can be viewed 
as the consequence of the Pauli exclusion principle and the flavour-exchange symmetry 
described above. Interactions favour such conhgurations in which different fermions do 
not overlap, and as a consequence particles are promoted to higher single-particle levels 
respecting the Pauli principle for fermions with the same spin. At the same time, the 
ground-state of the system has to remain in the symmetric conhguration under exchange 
of two particles with opposite spins (this is related to the antiferromagnetic structure of 
the ground state, see [47]). These two requirements force the system to excite favours 
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Figure 3. Single-particle densities (thin blue line, heavier flavour) and Pi{x) 

(thick red line, lighter flavour) calculated in the ground-state of the system for different 
numbers of fermions with a = 40/6 and strong interaction g — 10. In contrast to the 
previous results obtained for a = 1, in this case the separation always occurs in the 
lighter species, independently on the distribution of fermions between flavours. In 
particular, the separation is present also for an equal number of fermions = iV^. 


with larger iVo- to higher states than the particles of the minority. In consequence, 
the probability of finding a minority fermion in the centre of the trap decreases slower 
than for majority fermions. This is the reason why this kind of separation requires an 
imbalanced particle number. It is worth noticing that this argumentation is valid only 
for the ground-state of the system which in the limit of inhnite interactions becomes 
degenerate with states of other exchange-symmetries. Those other states may not 
manifest separation in any of the flavours [47]. 

For a 7 ^ 1, the separation has a different origin. In this case, the single¬ 
particle states (f)ia are not symmetric under flavour-exchange and as a consequence, 
the interaction energies between the two particles occupying different levels of the 
harmonic oscillator are not symmetric either. This fact leads directly to the separation 
in the flavour with smaller mass. To explain this mechanism, let us concentrate on 
the case = 3 and = 1. In the limit of vanishing interactions, the ground- 











































state of the system has the form |210;0) = 1/a/6 a^^aj^^ag^aQ^lvac). When interactions 
are present then particles of different spins try to avoid each other to minimize the 
interaction energy and as a consequence they change their Fock-space conhguration. 
The interaction energy can be minimized by exciting particles to higher single-particle 
states. In the model studied here, the energy cost of such an excitation does not depend 
on the flavour. However, due to the violation of the flavour-exchange symmetry in the 
interaction sector, the energy gain is larger when lighter particles (i.e. particles with a 
larger spread of the wavefunction) are promoted. In consequence, the Fock states of the 
form |210;k) = l/\/6 a 2 .|-ci|.|-aot®l 4 ,l^^^) with fc > 0 start to contribute and dominate 
the exact many-body ground-state of the system. Direct inspection of the ground-state 
obtained from the exact diagonalisation method proves indeed that this simple heuristic 
argumentation is correct. For example, for g = IQ and a = 40/6 the projection of the 
true ground state |Go) onto the non-interacting ground-state is |(210;0|Go)P < 10““^. 
At the same time, the cumulative probability of hnding all three heavy particles in the 
lowest states, i.e. I(210 ;/c|Go)P ~ 0.9. This means that the energy of the ground- 
state is minimised mainly due to the mobility of the lighter particle, i.e. the particle 
with the larger spreading of single-particle states. 

We have checked that the mechanism described above is very general and does not 
depend on the number of fermions of different flavours. In Fig. 3 we present the single¬ 
particle densities for different number of fermions and obtained in the strong 
interaction regime g = 10. In all plots, the mass ratio is a = 40/6. It is seen that the 
separation induced by a mass difference always occurs in the lighter component (red 
thick line). Our results performed for all conhgurations with + < 10 conhrm that 

the lighter component is always divided into two domains regardless of the number of 
fermions. 

It is worth noticing that for = 1 and = 2 our results obtained in the 
strong interaction limit are consistent with the density prohles obtained in [39] where a 
variational approach was used for this three-body problem. 

6. Non-zero temperatures 

The method used allows us to study the properties of the system also at non-zero 
temperatures, i.e. in a mixed state of the form px = >2“^ ^•pi|Gj)(Gi|, where pi = 
exp(—£’j//cs7") and Z = YhiVi is the partition function. A discussion of whether the 
one-dimensional system described by the Hamiltonian (1) reaches the thermal state px 
when interacting with the thermal bath is beyond the scope of this article. However, 
our aim is to discuss the influence of the interaction with an environment on the 
phenomenon of phase separation. From this point of view, the temperature T can 
be viewed as an effective parameter which controls strength of the interaction with an 
external environment. 

After dehning the mixed state px, the whole analysis can be performed in analogy to 
the zero temperature case provided that the single-particle densities Pa{x) are calculated 
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Figure 4. Single-particle densities (thin blue line, heavier flavour) and Pi{x) 

(thick red line, lighter flavour) calculated in the thermal state pj with = Z and 
= 1 for g = 10 as a function of temperature T. For a = 1 the phase separation 
induced by the imbalanced number of fermions is instantly destroyed at non-zero 
temperature. In contrast, for a = 40/6, when the phase separation is driven by a 
mass difference between fermions, the single-particle density of lighter flavour is almost 
unchanged and it is very robust to the mixing with higher many-body states. 


in the state pi, i.e. 

p„(x, T) = Tr p^ '!• 


( 6 ) 


In practice, we cut the number of exact many-body eigenstates at the number M (dehned 
in Sec. 3) which is temperature dependent. For a given temperature, this number is 
determined from the condition that the hnal single-particle densities do not change 
signihcantly when M is increased. It is natural that in the limit T —)■ 0 we reproduce 
the previous results obtained for Po-{x) dehned by eq. (5). 

In Fig. 4 we visualize the inhuence of non-zero temperatures on the phase separation 
for a number of particles equal to = 3, = 1 and strong interactions g = 10. As seen 

in the upper panel for a = 1, the phase separation is destroyed almost immediately when 
the higher eigenstates of the Hamiltonian contribute. Density prohles are becoming 
smoother and spatial density huctuations vanish. In contrast, for a = 40/6, the 
separation induced by a large mass imbalance is almost insensitive to the temperature, 
i.e. the lighter particle remains in the Schrodinger cat-like state. The only visible 
consequence of higher temperature is a smoothing of the bulk Favour density. This result 
confirms that the mechanism of the separation induced by a mass imbalance between 
Favours has a completely diFerent origin than the separation induced by interactions 
for a system with Favour-exchange symmetry. 


7. Conclusions 

We have studied the properties of a two-component mixture of fermions in a one¬ 
dimensional harmonic trap for diFerent mass ratios between atoms. We show that the 
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mass difference induces a phase separation in the lighter flavour independently of the 
distribution of fermions between the flavours. The mechanism observed is general and 
is forced by the interaction energy cost when particles are excited to higher states of the 
harmonic oscillator. Finally, we show that in contrast to the standard separation induced 
by an imbalanced number of fermions, the separation driven by the mass difference is 
robust to non-zero temperatures. We believe that our predictions can be conhrmed in 
upcoming experiments with ultra-cold mixtures of fermions. 
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